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Abstract 

We propose a model of Dirac neutrino masses generated at one-loop level. The origin of this 
mass is induced from Peccei-Quinn symmetry breaking which was proposed to solve the so-called 
strong CP problem in QCD, therefore, the neutrino mass is connected with the QCD scale, Aqcd- 
We also study the parameter space of this model confronting with neutrino oscillation data and 
leptonic rare decays. The phenomenological implications to leptonic flavor physics such as the 
electromagnetic moment of charged leptons and neutrinos are studied. Axion as the dark matter 
candidate is one of the byproduct in our scenario. Di-photon and Z-photon decay channels in the 
LHC Higgs search are investigated, we show that the effects of singly charged singlet scalar can be 
distinguished from the general two Higgs doublet model. 
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I. INTRODUCTION 



Small quantities arising in physics usually requires the use of new symmetries for ex- 
planations A good example of these is the Peccei-Quinn (PQ) symmetry that plays 
the role for the solution of the strong CP-problem j^, 3|, in which a ^-arigle which appears 

2 

in the QCD lagrangian, Cqcd G ^3^^%^"'^'^^ CP is violated [4-|6|. The instanton 
solution to the gluon field equations satisfies n = J d'^xG^^^G"'^^ with n's are integers 
and represent topological charges Ql- The QCD vacuum state hence can be parametrized 
as \6) = Yln=°^oo^^"^^\''^) ^ periodic with period 27i. Furthermore for nonzero quark 
masses the chiral anomaly relates the weak phase in quark masses to the QCD 6'-term, one 
can parametrize the 6 angle as ^ = 6 — arg{detmq} [sljO]. It induces a neutron electric dipole 
moment (EDM) 1CH14| and the current experimental upper bound set the best constraint 



on 6 to be smaller than 0.6 x 10" 
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15| . This extremely suppressed quantity is called the 



strong CP problem. The Peccei-Quinn solution to the strong CP problem postulates a global 
chiral U{1)pq symmetry and makes 6 a dynamical variable, and the shift symmetry of the 
Nambu-Goldstone boson, axion, corresponding to U{1)pq Q, will set 6 zero at classical 
potential 18[ . At one-loop level the chiral anomaly will break the shift symmetry and as a 
result the axion is not massless but requires a small mass 



~ {mu+nid) fa 



A2 
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Here is the U{ 
and outer space 



)p n b reaking scale and f^^ is the pion decay constant. The laboratory |22| 



23 



27| searches have set lO'^ GeV < /a < 3 x 10" GeV as the allowed 
regions, therefore, the axion window is 3 x 10~^ eV > nia > 10~^ eV. 

On the other hand, another small quantity that puzzles high e nerg y physicists is the 



masses of neutrinos measured from neutrino oscillation experiments |22|. The key point to 
understand neutrino physics lies on whether the neutrinos are Dirac fermions or Majorana 
fermions. This ambiguity comes from the fact of zero electric charge carried by neutrinos. 
The tiny neutrino masses may be explained in terms of lepton number (L) symmetry which 
is a global U{1) quantum number tagged on lepton sectors in the standard model (SM). 
If U{1)l is broken one can write the dimension-5 Weinberg operator to generate neutrino 



masses m^, oc 



HHLL 



here H and L are the SM Higgs and the left-handed lepton fields 



respectively, and A^ is the breaking scale of U{1)l. In this case neutrinos are regarded as 



1 The are the SU{?,)c gauge fields with a = 1, 2, 8 and G"^'" = \e^"'°'f^Gl^. 
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Majorana fermions. However, the Majorana or Dirac nature of the neutrinos is unknown and 
is awaiting for the experimental determination from some lepton number violating processes 
such as the neutrinoless double beta decay. It is important to consider the possibility that 
Dirac neutrino masses may also connect with some global symmetry and how those small 
quantities we observe in physics are related to each other. In this paper we propose a simple 
Dirac neutrino mass model which is generated by PQ symmetry breaking, and hence the 
neutrino masses are closely related with axion mass. 

This paper is organized as follows : in section II we propose the Dirac mass model which 
is embedded with PQ symmetry. Section III we consider leptonic rare decays and neutrino 
oscillations data to investigate the parameter space of the model. Some phenomenological 
implications to LHC Higgs search, dark matter, and leptonic flavor physics of this model 
are discussed in section IV. Then we conclude our results in section V. 

II. THE MODEL 

Particle content and their quantum numbers of the model are listed in Table [B The 
Y, L and PQ represent the hypercharge, lepton number, and U{1)pq charge respectively. 
Two Higgs doublets are introduced since the existence of the PQ-symmetry one need two 
independent chiral transformation for the up-type and down-type fermions. We consider 
the scenario of neutrino is Dirac fermion three right-handed neutrinos i^/j. (i = 1 — 3) 
are assigned to our model, and hence the theory is lepton number conserved, sf and 
are SU (2) ^-singlet charged scalars, and a is the axion field which is the Nambu-Goldstone 
boson of the spontaneously broken U{1)pq. Notice that z/^j's are complete neutral under 
gauge symmetries and PQ symmetry and only carry the L quantum number. The Dirac 
neutrino mass term is forbidden by the PQ-symmetry at the tree level and is generated at 
one- loop level after the PQ symmetry breaking by utilizing the charged scalars sf and s^. 

The new Yukawa interactions of the model for leptons are given by 

C = y^p{LL^)lR^Hi + f^pLlJa2{LL^)st + KiI'r^vrA + h.c, (1) 

where c denotes charged conjugation, = e,fi,T, and cTj (i = 1 — 3) are the Pauli 
matrices. In general, y and h are complex matrices, and / is an antisymmetric matrix due 
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TABLE I: Quantum numbers of U{1)pq and gauge symmetries for leptons and scalars. 
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FIG. 1: Induced Dirac neutrino mass 

to the Fermi statistics. One can choose the basis of leptonic mass eigenstates Ll, such 
that ?/ is a diagonahzed matrix. Also fij can be chosen to be real by rephasing and by 
transferring the phases into Ir, therefore only hij are complex. The scalar potential can be 
written as 

+ \i{H\Hi){HIHi) + \2{HIH2){HIH2) + \^{H\Hi){HIH2) + HH\H2){HIHi) 

+ {H\Hi)[h\si\'' + + /a|a|2)] + {HlH2)[gi\si\'' + g2\s2\^ + ga\a\''] 

+ hi\si\'^ + h2\s2\^ + ha\af + /l3|SinS2p + |ap(/lai|si|^ + ha2\s2\^) 

+ [K{s'{s'^f + h^{HlHi)o? + ^s^s'^ a + h..c.]. (2) 

All parameters in the potential are real though /14, and /i are in general complex param- 
eters, one can absorb their phases by redefining (s|s2), (hIHi) and a respectively. For the 
invisible axion, f„, (7^, hai, ha2, ^5, should be very small to make other scalar mass not 



too heavy [29l432l| and are irrelevant in our discussions. 

The leading contribution to the Dirac neutrino mass at one-loop level is shown in Fig. [1] 
and is given by 
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{M^)ai = - p faphpil{ml^,ml^, ml) , (3) 
oTT mp 



where /(m^^^, m^2, is defined as 

Tf 2 2 2\ '^l t T^'il 1 'mil 1 "^?2t /.\ 

/ 771^^,771^2, m| = ^ o o ^Qg^ 2 2 ^Og ^ " 4) 

777fi - 777^2 mj^ - mj mj 777^2 " mj mj 

In the hmit of msi^2 ^ "7^, the neutrino mass matrix is proportional to charged lepton mass, 
as given by 



{My)o,i = ---^fal3hi3imi3^ ^-^log — 

OTT^ 777;i — 777^0 777 



si 



^s2 "^s2 



1 , , f^mp Aqcd , mil nf u fK\ 
^Japhpi 2" log = -Cfapmphpi. (5) 



OTT^ 777a "^sl ~ ''^s2 "^s2 

We have replaced the PQ-symmetry breaking scale by the axion mass and the and the 



QCD scale in Eq. (|5]). Note that there are similar construction discussed in Refs. [33|, l34j . 
however, our setup shows the interesting relation between the neutrino Dirac mass and QCD 
dynamics. Furthermore, the axions can be the dark matter candidate constituting of 25% 
of energy density in our universe. We can see that the large value of fa lifts up the mass 
scale of Dirac neutrino masses and in order to keep the smallness of Mj^ a combination of 
suppressed factors will be needed such as the loop factor, the Yukawa couplings /, /i, the 
charged lepton chirality suppression, and the parameter /i. Now we roughly estimate the 
scale of /i which is the key parameter controlling the overall mass scale of Dirac neutrino 
masses and a investigation of parameter space will be discussed in section III. For M^, ^0.1 
eV, 777,, 2 ~ C(100 - 1000) GeV, and fa ~ 10^° GeV, we have / ~ lO'^, h ~ lO'^, and 
/i ~ 0(100) keV. Let's make two comments on the low scale /i and Dirac neutrino masses : 



1. We can explain the small /i by implementing the Froggatt-Nielsen mechanism 35| with 
PQ-symmetry. For example, if we assign the PQ quantum number of a as |- in Table [U all 
the terms in the Lagrangian will not change except the /i-coupling in the potential. One can 
write the effective operator as 

1 (M)-2(a)s+s2-a, (6) 



A""2 U(l)pQ breaking fa 

thus /i = (^)"-2(a) at low energy scale and can be tuned to a small quantity. Here 
the PQ symmetry can be broken dynamically by some condensates of a new technicolor- 



like interaction at high scale 



36 



37|. 2. The mass scale of Dirac neutrino masses is not 



necessarily small if neutrinos are Majorana fermions. Although we consider neutrinos are 
the Dirac fermions, the main concern in this paper is the origin of the Dirac neutrino 
mass, which in general does not forbidden the possibility that neutrinos are the Majorana 
fermions. Therefore one can still have heavy right-handed Majorana masses as inspired by 
the grand unification theories and obtain small neutrino masses through canonical seesaw 
mechanism. The goal in this paper is to point out the Dirac neutrino mass is generated by 
the PQ-symmetry breaking. 

III. CONFRONTING WITH NEUTRINO OSCILLATION DATA AND LEP- 
TONIC RARE DECAYS 

From the standard formalism can be diagonalized by 

= VlunsMuV^R^ = -V^MNsiCfMi,^^^h)V;^\ (7) 

here Vpmns is a unitary 3x3 Pontecorvo-Maki-Nakagava-Sakata (PMNS) matrix and Vr 
is the transformation matrix for right-handed neutrinos. For convenience we define FH = 
^PMNS^i/d^ag "with H = Mi^^^^{—hV^) and F = Cf. Due to the anti-symmetric nature of / 
matrix, the lightest neutrino is exactly massless in this model. Therefore, we can multiply 
a transformation matrix A to both sides of the mass matrix to reduce one row in the left 
hand side of Eq. ([7]). One obtain 

rH = AVpMmM.,^^^, (8) 

where 
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Then we have 



Per H el + F^rH^i FeT-ffe2 + F^rH^2 FerHeZ + Ff^rH^s 


\ 



"^i^lV^l "1^2V"^2 "Ij^gV^a 

-m^lKl -"^;.2K2 -"2^.3^3 

with Fi = (F^,,-Fe,,Fe^). 

For the normal hierarchical spectrum with m^.^ = we have the following relations : 



(10) 
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The other terms give H^2i H^^^, Ht-2, H^-^ in terms of F^,-, Hf,2 and i^es 
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(13) 



The real F matrix leads to the dependence of two phases in Vpmns and only one irreducible 
phase remains. 

Similarly for the inverted hierarchical spectrum, m^,^ = 0, we obtain 



F = — F F 
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-Vr2 , Ht-2 = rr—He2 + 
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/i2 , 



eS 



MT 
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U TJ t3 TJ 

J^e3 , -HrS — -TT-rt, 
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(14) 



eS 



We will use the central values of the most recently global fitting of the neutrino oscillation 
measurements 38|] (see Table [Tll) in our analyses. In the meantime the appearance of the 
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sin^ 6112 sin^ 6^23 


sin^ 6*13 






0.30 ± 0.013 0.4llJ]:[!^^ 


0.023 ± 0.0023 


(7.50±0.185)xl0~5eV2 


(2.47t°;°^?) X 10-3eV2 



TABLE II: Neutrino oscillation data. 



new scalars will provide the extra contributions to the lepton flavor violation processes. We 
investigate the parameter space of the model in terms of the constraints from leptonic rare 
decays in the following. Before that it is worth mentioning that the size of parameter C is 
proportional to the factor fifa/{m1i — m^2) log('^si/''^s2)- If there exists a large hierarchy 
between m^i and 771^2, say, 771^2 = krUsi with k being a large factor, C is approximately 
inverse proportional to m^2- However, the positive mass eigenvalue condition msims2 > P-fa 
will also lead to the result that the upper bound of C is proportional to 1/k. In general C 
can not be too small in order to keep the perturbativity of the Yukawa coupling / and h. 
So the largest allowed hierarchy between m^^ and is around k ~ O(IO^). We will also 
discuss the implication to the parameter space with a hierarchy scenario. 



A. /i — )• 67 

The two relevant Yukawa interactions providing the flavor violations in charged lepton 
sector are given by 

jC = -2fapll^l'LfiSl - hail'j^J^R,S2 + h.C. (15) 

Without loss of generality here z/r. is the mass eigenstate and we absorb the mixing matrix 
— Vr into the coupling h. In SM + z/^j the one-loop contribution is constrained stringently by 
the Glashow-Iliopoulos-Maiani (GIM) mechanism. In this model the main contribution is 
the one-loop diagrams with photon emission attached to the charged scalars s\ or s\ in the 
loop. Currently the latest result from MEG collaboration gives — 67) < 2.4x 10~^^ 39 |. 
For la ^ 1 131 the effective Lagrangian can be generally written in the form 

C = -hpCT,,iAnPR + AlPl%F^\ (16) 
where Al,r for this model is given by 
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Normal Hierarchy 



c=io-: ■ 

C=10"= 
C=0.9x10"* 




Inverted Hierarchy 




FIG. 2: The allowed region in h^i — f^^ plane in (a) the normal hierarchy, and (b) the inverted 
hierarchy. Different curves are generated with different C (blue: C = 10~^, green: C = 10~^, red: 
C = 0.9 X 10~^). The other parameters are taken as nisi = = 500GeV, /ie2 = /^es = 0.5. 



Note that the hmit ma ^ 2 have been apphed to the above formula. The decay rate of 
la — ^ Ipl is — )■ /^7) = (m^/167r)(l — m^/m^)'^(|y4ijp + One can compare it with 

the lepton three body decay width T{la — )■ If^Va^i^p) = G^pm^/ 19271"^. The branching ratio 
/i — 67 in our model is obtained by 



B(/i ^ 67) 



T{fi 67) 



Of. 



7687TGI 



16E,l/e7lVM7P , E^\he.\'K.\ 



m 



si 



(18) 



The results are given in Fig. |2l Here we consider f^^ is positive for the normal hierarchy 
spectrum and f^r is negative for the inverted hierarchy spectrum, respectively. The allowed 
regions are rather small and sensitive to the value of C . For example, the parameter space 
will shrink to zero if we take C = 0.8 x 10~^ with the same inputs. 



B. ^ — e conversion 



The one loop diagrams to fi — e conversion include photon penguin diagrams, Z penguin 
diagrams, and box diagrams. Again the contribution from the diagrams involving W boson 
exchange is suppressed by GIM mechanism and the leading contributions come from the 
penguin diagrams with the charged scalars sf and in the loop. In contrast, the leading 
Z penguin contribution is suppressed by the light charged scalars. Therefore only photon 
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penguin needs to be taken into account, the corresponding effective Lagrangian is given by 

V ^ si j s2 g 

In the above we used the shorthand notation sw = smOy/. The fi — e conversion with 
nucleon in atom has been calculated in details in Ref. 40|, and we will adopt their notation 
in what follows. The general interactions associated with fi — e conversion are written as 



C 



cff 



V2 



-[m^ea^'iARPn + ALPL)f^F^u + h.c] 



G 

■--^[e{gLsig)PR + 9RS{q)PL)Mq + e{gLP(q)PR + gRP{q)PL)lJ^qibq + h.c] 
G 

■^[K9LV{q)l^PL + 9RV{q)l^PR)Mli,q + e{gLA(q)l^PL + gRA{q)Y'PR)f^q7,i75q + h.c] 
G 



^ y^(i{9LT(q)(^'"'PR + 9RT(q)(y'"'PL)M(^tiyq + h.C 



(20) 



where Ar^l is related to the dipole interaction with photon, gs,p,v,A,T indicate scalar, pseu 
doscalar, vector, axial vector and tensor couplings, respectively. Comparing the above for- 
mula with Eq. ( |T9l) . we have 

_,(<?) _ Qq W^2 (y^4/*^/^a^) , guv = (^^ V ~ 



QC 2 "''W\ 
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(21) 
(22) 



1927,2^2 vz_. —^2^- 192vr2^2 _ 

and other couplings vanish. The rate of fi — e conversion with nucleons in atom A is usually 
normalized to the rate of muon capture by A. The conversion-to-capture ratio can be derived 
in the form 





2Glml 


capt 


with g^^Rv 


- 2(7^"^ 
~ '^yL,RV 



(leAnD + -g^lv^^^ + c,tiv^-^\' + \eA,D + ~gtlv^^^ + -g^^lv^ 



(23) 



fl'["ra^= &y + 25'l W- P' and V^^f'") are overlapped 



fi — e conversion in different nuclei 41 
C. fi — 3e 



functions which can be found in Ref. 40 1. The constraints from experimental results for 



45[ are weaker than what's given by /i — )■ ej. 



Penguin diagrams and Box diag rams contribute to this process. The experimental upper 
bounds is i?(/x — > eee) < 10^^^ [46|]. SM with right-handed neutrino singlets contribution to 
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this processes is suppressed by the neutrino mass. The corresponding one loop diagrams in 
this model for yU — )■ 3e are similar to those for fi — e conversion, with the quarks replaced by 
electrons. In general the effective Lagrangian for n — )■ eee is 



+eY{aLPL + aRPR)ee-ff,PLfi + eYibiPL + bRPn)ee-f^PRfi] . (24) 



£(/i eee) = ;=-[e7^e^ezo-'"'(87rae"i^)(^i?-Pi? + ^l^l)/^ 



The branching ratio can be calculated as 

B{i2 eee) - (a^ + bl) + 2{al + b^) - 32iTae{AR{aR + 2aL) + AL{bL + 2fe/j)) 

+256n'ali\AR\' + \A,\')(^A\n^ - ^) , (25) 

where the parameters ql^r and bi^R are given by 

1 1 

aL,R = -zr^mlriJ2^f*eaUc.-^) ^ud = - T-^m^ (^ Z^:, V (26) 



from the photon penguin diagrams. While for the box diagrams the leading order of a]^ 



box 



and 6^™ vanish, we have 



box ''^w 



E4/M-/eJlE4/ea'/e.'| (27) 



and 



^^""^ ~ Q9^2Xr,2 I E V^eill E ^ei'^^i'l (28) 

respectively. As shown in Fig. [2] the Yukawa coupling / is in the range of (9(10~^~~^), we 
can safely ignore the box diagram contributions in — )■ 3e decay. Therefore, the parameter 
space is looser to those of /i — e7. 

In the limit of mass hierarchy between rrisi and 772^2, that is, m^i > 777,^2 or 777^2 > ^si 
cases. The parameter space for Yukawa couplings h and / in both normal hierarchy and 
inverted hierarchy neutrino mass spectrum respectively are shown in Fig. [3] and Fig. HJ Here 
we take rrisi = 20777^2 and 771^2 = STTi^i as the reference points. In these cases, we are only 
able to give a severe constraint on one of the couplings, h ot f, and we illustrate our results 
by taking 500 GeV mass scale to the lighter scalar field. 
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Normal Hierarchy, mj,j>m^T Inverted Hierarchy, m^i>m^2 




FIG. 3: The allowed region in h^i — f^r plane for the case of m^i = 20ms2 in normal hierarchy and 
inverted hierarchy neutrino mass spectrum respectively. Blue line corresponds to C = 10~^ and 
Red line represents C = 0.9 x 10~^. 



Normal Hierarchy, m^2>"r^i Inverted Hierarchy, m^2>"isl 




FIG. 4: The allowed region in h^i — f^r plane for the case of 771^2 = Sm^i in normal hierarchy and 
inverted hierarchy neutrino mass spectrum respectively. Blue line corresponds to C = 10^'^ and 
Red line represents C = 0.9 x 10~^. 

IV. PHENOMENOLOGY 



A. electromagnetic moments of leptons 



1. muon g — 2 



The anomalous magnetic moment o 
SM. The current experimental results 



the muon is an observable as a precision test to the 



47 



481 ] reported that the muon g — 2 deviation from 
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the SM prediction is a^(exp) — a^(SM) = (287 ±63 ±49) x 10~^^, which indicates some new 
physics contribution. Since our scenario is essentially a two Higgs doublet model (2HDM) 
plus two singly charged scalars sf and . The anomalous g — 2 oi muon in this model is 
given by 



a 



s2 



.K< , 22tan2/3 tan^ /3 

COS"' p COS"' p 

where and if are scalar particles with h being the SM-like Higgs, as well as the pseudo- 
scalar field A and the charged scalars . a, /3 are the mixing angles (s/3_q, = sin (/3 — a) 
and c^_Q, = cos (/3 — a)) defined in the usual 2HDM (see e.g. [49| and references within). 
Comparing with the current experimental result, our model gives a rather small contribution 
to Aa^ but still within the experimental errors. 



2. magnetic moments of neutrinos 



Neutrinos can have magnetic moment when they are massive. The present upper bound 
of neutrino magnetic moment from experiments is /i^ < 3.2 x 10 fiB |[50]. The ordinary 
leading order contribution in SM including ri ght- handed Dirac neutrino is the exchange of 
W, which is around (3.2 x 10~^^yUB)(mjy-/eV) [5ll] with the Bohr magneton fiB = e/2me. In 
this model, the main new contribution comes from the mixing of loop since the 

same diagram also generates neutrino masses, given by 



l/i^l = 2mie(log 



m 



s2^ 



mz 



) ^ IO-^Vb- (30) 

Note that it only depends on the neutrino and scalar masses, and this contribution is com- 
parable with that associated with W exchange. It is understandable that without imposing 
a symmetry [52| or employing a spin suppression mechanism to rrii, 53|, the generic size of 
the Dirac neutrino magnetic moment can not excess 10~^^/iB M]- The contributions from 
other charged scalars such as H~^, sf, and in the loop are less than O{10~'^^)fiB as a 
result of they do not contribute to neutrino masses. 
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3. EDM of leptons 



Since there is a physical CP phase among the Yukawa couphngs / and h, in general, we 
have the new contributions to the electric dipole moment of charged leptons. The exchange 
of singly charged scalars gives the EDM to the charged leptons at one-loop level as 

= -7AEM2/^r.(/nO^^[— 5 (^) - —5 (^)] , (31) 
(47r)^ ^ mj.^- 17112 i^si V"^si/ \m1.2j 

with the function B{x) is defined by 

2 

Therefore we conclude that the extra contribution is smaller than 0( "^"2 )1 which means 
Mel ^ 10~^^ecm. Similarly, the neutrino EDM generated by mixing in the loop gives 
\dj\ < 10~^^ecm. Though nonzero but both are unobservably small. 



B. dark matter 



It is known that the axion field can be a dark matter candidate constituting a significant 
fraction of energy density in our universe. For the purpose of completion we briefly review 
some aspects of this scenario in this subsection. The properties of the invisible axion are 
determined by the breaking scale fa where its mass and the interactions are inverse pro- 
portional to fa- Hence the invisible axion is a very light, very weakly interacted and very 
long-lived particle. Axions with mass in the range of 10~^ — 10~^ eV were produced during 
the QCD phase transition with the average momentum of order the Hubble expansion rate 
(~ 3 X 10~^ eV) at this epoch and hence are cold dark matter (CDM). Their number density 
is provided by 

1,0.6 X 10"^ eV, 7,0.7,. 

h is the current Hubble expansion rate in units of lOOkms^^Mpc"^. Here we assume the 
ratio of the axion number density to the entropy density is constant since produced and 
the contribution from topological defects decay is negligible. Recently it was pointed out 
that the CDM axions would form a Bose-Einstein condensate due to their gravitational 
interactions 55|. Furthermore the rethermalization process is so fast that the lowest energy 
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state of the degenerate axion gas consistent with a nonzero angular momentum, as a result 



that a "caustic ring" structure may form in the inner galactic halo 



56- 



The feature 



would make the axion a different dark matter from other CDM candidates, we refer readers 



to the references 



55 
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for details. 



C. h 



77 



We closing our discussion on phenomenology with investigating the LHC Higgs results. 



Both ATLAS 



59| and CMS 60| have announced the discovery of a new boson at a mass of 



125 GeV which is consistent with the SM Higgs boson via the combined analyses of the 77 
and ZZ channels. However, the precise values of both production cross sections and decay 
branch ratios of the new resonance need to be measured to compare with those predictions 
from the SM. It was pointed out that the branching ratio of Higgs decay into two photons 
excess about 1.56 ± 0.43 and 1.9 ± 0.5 times than the SM prediction in both CMS 



ATLAS 



60| and 



59| collaborations data. Although the deviation is still within the SM expectations 



at 2a level, one may consider whether there are new physics effects (see e.g. 6l|-[70| and 
references therein), in particular, new scalars particles have been widely treated as possible 
sources. We study the implications of Higgs to di-phonton and Higgs to Z-photon decay 
channels in our model. For simplicity we just consider one singly charged scalar and omit 
the subscript in this subsection. The SM Higgs production cross section is modified by the 
additional doublet scalar in our scenario, 



-(^^^ + Sa)Ai/2{Tt) - -{CatanP - Sa)Ai/2{n) 



(34) 



128^271 

notice that the bottom quark contribution is not negligible due the enhancement of large 
tan/3 and Ai/2(r) = 2[r + (r - l)/(r)]r-2 with n 

f{r) 



(sin-' v^)2 

-i[iogi±^i^ 



The function /(r) is defined by 
r < 1 

(35) 



,r< 1 
? ,r>l 



In type-II 2HDM the decay rate of /i — )■ 77 is given by 

A' A 

+s^Ai{tw) + -^Ao{th+) + -Ao(r,+ )|^ 



(36) 



15 



where Ai(r) = -pr^ + 3r + 3(2r - l)/(r)]r-^ and Ao(r) = -[r - /(r)]r-2. The coef- 
ficients A and A' are defined as A = vfi^s+s- /nT-l and A' = vfinn+H- /^"H+y where ^hs+s- 
and Hhh+h- are the related trihnear couphngs in the Higgs potential One can see that 
the effects of doublet scalar charged component and the singly charged singlet are 
indistinguishable. The reason is the lack of knowledge of scalar potential. For h — j- ^7 
channel the decay width is written as 

^^7 - \^-t^zMh) X 

+ ^^v4i/2(r/t, ftt) + (c„ tan/3 - ^^^1/2(77;,, k^) 



tanp cvi/ cvK 

1 ' 
+s„v4i(?7h/, kw) + (o - ) A' v4o - s^Av4o(?7s+, 



(37) 



.2 „, ,2 



where Ai/2(r/, k) = /i(r/, k) - 12(77, k), Ai(r/, k) = Ciy{4(3 - 'j-)h{j]. «:) + [(1 + f )^ - (5 + 
^)]/i(?7, k)}, and Ao{ri, k) = h{ri, k) with the parameters rfi = and Ki = Functions 
/i(?7, k) and hiv^ are given as 

and 

^2(r/, /.) = - /('^"')] (39) 

with 



9W = I (40) 

respectively. Again we see from the Eq. ( l37l) that the singly charged singlet effect is hidden 
in the 2HDM. In order to extract the information of from the decays we found the ratio 
of Aq{t], k) /Aq{t) lies in the range of 0.50 — 0.53 with charged scalar masses above 100 GeV. 
We then take the ratio as a constant and subtract the contributions in F^^ and F^^. 
The results of the effects in terms of the parameter cos a with tan /3 = 10, 50 are shown 
in Fig. [5l Here we take the values i?^-v = -2^2— = 1.5 and Rz^y = "^'^ = 1 with an 
and (Jiigj^i are the production cross sections for the Higgs to ii channel in our model and 
SM respectively. Since the new physics effects can have both constructive and destructive 
interferences with the SM ly^iy^ amplitude, there are four lines exist in Fig. [51 The 
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0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

cosa cosa 

FIG. 5: Solutions of s"*" effects in the terms of cosa by assuming AQ{rj, k)/Ao{t) = 0.5 with 
R^^ = =1.5 and Rz^ = "^^^ =1. an and aua,, are the production cross sections for the 

Higgs to ii channel in our model and SM respectively. 

results of effects with different values of Rz-^f and the ratio of v4o(?7, h)/Aq{t) for the fixed 
R^^ = 1.5 are shown in Fig. |6l Here Rz-y = 0.5, 1, 2 and Aq^t], k)/Ao{t) is taken to be 0.53 
and 0.50 respectively. It shows that taking the ratio of Aq^t], k,)/Ao{t) as a constant is a 
good approximation in the regions of parameters we are interested in and how the deviation 
of Tz-y from the SM prediction for the singly charged singlet scalar s~^. Finally, we estimate 
the background contribution for h — WW* — eufiu channel from s^s^ production. Both 
singly charged scalars are off-shell and we found the partial decay width for the center of 
mass in one of virtual scalar to be around rriw is negligible compared to the SM prediction. 

V. DISCUSSIONS AND CONCLUSIONS 

We investigate a model that neutrinos are Dirac fermions and their masses are generated 
from the Peccei-Quinn symmetry breaking at one-loop level. As a result the neutrino mass 
is related to Aqcd and axion is appeared to be a good candidate of dark matter in explaining 
the missing energy density of our universe. Leptonic rare decays constrain the model pa- 
rameters severely, therefore, the model can be tested in the near future. We also studied the 
implications of the new scalars to the electromagnetic moments of leptons and recent Higgs 
signals at LHC, specifically the /i — )■ 77 and h — Z7 decays. Finally we would like to make 
a brief comment on the Majorana extension of our scenario. Lepton number symmetry is 
one of the key to understand the underlying neutrino physics. The smoking gun signals to 
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< 

'< -10 























\Rz-ri 































FIG. 6: Effects of singly charged singlet scalar with values of Rz-y = 0.5, 1, 2 and Ao{r], k,)/Aq{t) = 
0.53,0.50 respectively for the fixed i?^^ = 1.5. Blues lines correspond to tan/3 = 10 and Red lines 
represent tan /3 = 50. Note that each figure corresponds to one of the lines exhibited in Fig. [5] 



resolve the question are OufSP decays and some lepton number violating processes at LHC 
Without the direct observations Majorana neutrinos and Dirac neutrinos are equally good in 
many aspects to describe the phenomena such as neutrino oscillations, leptogenesis, and Big 
Bang Nucleosynthesis, etc. Although we discuss the Dirac neutrino in this model and argue 
that the Dirac neutrino mass is originated from the Peccei-Quinn symmetry breaking which 
is well-motivated in QCD field theory. A Majorana masses of the right-handed neutrinos, in 
general, can be formed without violating any principle except the lepton number symmetry. 
The main modifications are that some terms to be included in the scalar potential such as 
H1H2S1 and HiH2as2 ■ Therefore, our discussions can be easily embedded the scenario of 
Majorana neutrino and diagonalize the neutrino mass matrix via the seesaw mechanism. 
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